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无 传感器 PMSM 中 基于 IGSO 优化 EKF 的 速度 估计 方法 
张 相 胜 ， 田 佳 文 ， 潘 “ 丰 


(江南 大 学 轻 工 过 程 先 进 控制 教育 部 重点 实验 室 ， 江苏 无 锡 214122) 


摘 要 : 为 了 提高 无 传感器 永 磁 同步 电机 (PMSMD) 控 制 系统 中 速度 控制 性 能 ,提出 一 种 基于 改进 群 搜索 优化 (GSO) 算 法 

pd 度 估计 方案 。 首 先 ， 分 析 了 PMSM 磁场 定向 控制 (FOC) 系 统 模型 ; 然后 ， 将 电机 的 d-g 轴 电 
、 电 流 和 转子 速度 作为 状态 变量 ， 构 建 EKF 中 的 状态 方程 来 估计 转速 和 负载 。 同 时 ， 为 了 提高 EKF 的 估计 性 能 ， 

We 通过 IGSO 算法 来 优化 EKF 中 的 噪声 协 方差 矩阵 QO 和 RR， 以 

此 获得 最 优 参数 。 仿 真 结 果 表明 ， 提 出 的 控制 系统 能 够 精确 估计 出 电机 转速 并 进行 有 效 控制 。 
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Method of EKF speed estimation for sensorless pmsm based on IGSO optimizing 
Zhang Xiangsheng, Tian Jiawen, Pan Feng 


(Key Laboratory of Advanced Process Control for Light Industry of Ministry of Education, Jiangnan University, Wuxi Jiangsu 
210094, China) 


IOS Abstract: In order to improve the speed control performance in sensorless permanent magnet synchronous motor (PMSM) 


control system, this paper proposed an extended Kalman filter (EKF) speed estimation scheme based on improved group search 
optimization (IGSO) algorithm. Firstly, it analyzed the PMSM field-oriented control (FOC) system model. Then, it used the 


motor's qd-d axis voltage, current and rotor Speed as state variables to construct the state equation, so as to use the EKF to estimate 


和 的- the speed and load. At the same time, in order to improvethe estimated performance of EKF, it estimated the square error integral 


(ISE) between the value and the actual value as the fitness function, and used the IGSO algorithm to optimize the noise 


covariance matrix O and R in EKF. The simulation results show that the proposed control System can accurately estimate the 
motor Speed and control effectively. 
Key Words: permanent magnet Synchronous motor; Speed estimation; extended kalman filter; noise covariance matrix; group 


search optimization 


reference adaptive system, MRAS)、 观 测 器 、 扩 展 卡 尔 曼 滤 波 器 

(extended Kalman filter, EKFP)、 模 糊 逻 辑 和 人 工 神 经 网 络 等 铝 。 
永 磁 同 步 电 机 (permanent magnet synchronous motor, PMSM) ”其 中 EKF 应 用 较为 广泛 ， 然 而 ，EKF 中 的 噪声 协 方差 矩阵 O 

有 具 有 功率 密度 高 、 效 率 高 、 鲁 棒 性 好 等 优点 ， 使 其 在 自动 化 、 和 R 的 取 值 直接 影响 估计 性 能 外 ,需要 取 最 优 的 参数 值 。 

交通 、 航 空 航天 等 领域 得 到 广泛 应 用 趾 。 对 于 PMSM 控制 ， 主 为 此 ， 本 文 提出 一 种 基于 改进 群 搜索 优化 (improved group 

要 分 为 磁场 定向 控制 (FOC) 和 直接 转 矩 控制 (DTC) 外 。 其 中 ， 都 search optimization, IGSO) 算 法 的 EKF 速度 估计 方案 , 并 通过 仿 


lf 
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是 通过 前 级 的 速度 控制 器 输出 信号 来 调节 电流 和 磁 通 。 为 此 ， 真实 验证 明了 提出 方法 的 有 效 性 。 提 出 方法 的 主要 创新 点 如 下 : 
速度 控制 成 为 整个 闭环 控制 的 关键 。 而 速度 控制 的 基础 是 准确 a) 将 电机 的 4! -9 轴 电 压 、 电 流 和 转子 速度 作为 状态 变量 ， 
获得 电机 的 实际 转速 器 。 构建 EKF 中 的 状态 方程 。 

传统 电机 控制 需要 速度 传感器 或 光学 编码 器 来 测量 转子 速 b) 对 传统 GSO 算法 进行 改进 ， 引 入 粒子 群 优化 (particle 


度 。 目前 , 为 了 降低 成 本 和 提高 PMSM 对 器 件 的 鲁 棒 性 , 通常 。 “swarm optimization, PSO) 算 法 中 的 全 局 最 优 和 个 体 最 优 融 合 技 
日 无 速度 传感器 的 PMSM 控制 外 。 其 中 ， 电 机 转速 是 通过 佑 ” ” 术 ， 改 进 GSO 中 搜索 者 的 位 置 更 新 机 制 ， 以 此 提高 其 收敛 速度 
计 模 型 来 预测 ， 估 计 模 型 通常 有 模型 参考 自 适 应 系统 (model 


并 
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1 PMSM 数学 模型 

PMSM 具有 不 同 的 4 轴 和 gq 轴 电 感 ， 这 取决 于 有 效 气 隙 的 
长 度 。 转 子 中 具有 极 的 部 分 由 4g 轴 定 义 ， 且 其 为 励磁 绕组 的 磁 
通 产 生 方 向 四。g 轴 定 义 为 与 4 轴 成 90° 电 角度 的 方向 。 图 1 所 
示 为 PMSM 的 结构 。 


图 1 PMSM 的 结构 


PMSM 通常 等 效 成 一 个 转子 旋转 坐标 系 (4 - 4 轴 ) 电 路 ， 如 


图 2 所 示 。 其 中 ，4 轴 和 4 轴 的 定子 电压 表示 为 
VW =ri, 4 -0@1 (1) 
di 
v=ri Pe -OH CO) 
gq 2 dt gq x 


其 中 : vy,、 ww、2i 记分 别 为 4-9 轴 电压 和 电流 ; 而 @,、x 
分 别 为 电机 转速 和 定子 电阻 。 和 和 分 别 为 4- 7 轴 的 磁 链 ， 
可 以 用 定子 电流 、 电 感 和 恒定 磁 链 4, 来 表示 : 
4 =Li th G3) 
4 =Li, (4) 
其 中 : L、 工 分别 为 4 -gq 轴 电感 。 将 式 (3)(4) 代 入 式 (1)(2) 中 
定子 电压 方程 可 修改 为 


l d ， 
Vi=rit+L, 二 — ,Li (5) 
pd Li 1 (6) 
Wy 二 gt a 十 ON 
电磁 转 矩 为 
3P ，， 要 
T 全 +(L, —L)ii,] (7) 
运动 方程 可 以 表示 为 
je (8) 
dt 


尼 系 数 ，0 为 转子 位 置 ，J 为 惯性 矩 。 


P 表 示 极 数 ， B 为 阻 


图 2 PMSM 的 等 效 电路 
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2 基于 EKF 的 PMSM 无 传感器 控制 模型 


2.1 PMSM 无 传感器 控制 模型 

对 于 PMSM， 本 文采 用 常用 的 磁场 定向 控制 (FOC) 系 统 ， 
将 PMSM 转换 为 单独 激励 的 等 效 直流 电机 。 图 3 为 PMSM 的 
磁场 定向 控制 系统 框图 ， 其 包含 了 脉 宽 调 制 PWMD) 道 变 器 、 
PWM 发 生 器 、PI 速度 控制 器 、 电 流 控制 器 ， 和 EKF 速度 /位 置 
估计 器 四。 


图 3 PMSM 无 传感器 FOC 控制 模型 


2.2 扩展 卡尔 曼 滤 波 (EKF ) 

卡尔 曼 滤 波 可 以 对 线性 系统 的 状态 进行 估计 ， 但 如 果 将 其 
应 用 于 非 线 性 系统 ， 则 会 非常 耗 时 站 。 这 是 因为 在 非 线 性 系统 
中 ， 状 态 变 量 的 方程 在 每 一 步 中 都 会 发 生变 化 ， 因 此 在 迭代 过 
程 中 不 能 进行 预计 算 。 这 一 不 足 可 以 通过 EKF 方法 来 解决 9。 
在 PMSM 数学 模型 中 ， 将 当前 y -g 轴 电 压 、 电 流 和 转子 速度 
选 为 EKF 的 状态 变量 , 那么 输入 变量 w 和 输出 变量 y 可 如 下 定 


义 : 

i . 
必 四 站 

O， _ 

寻 此 ， 离 散 时 间 系 统 方程 可 表示 为 

X= fr (Ke 1 We) (10) 
yx = (Kv) (11) 
其 中 :，x 为 状态 向 量 ， 为 输入 向 量 ，w 为 随机 状态 噪声 ， 


y 为 噪声 观测 或 待 测 变量 ，w 为 测量 噪声 。 

2.3 ”EKF 算法 的 步骤 描述 如 下 : 
a) 初 始 化 状态 向 量 x(0), P(0) 和 协 方差 矩阵 ， 即 O,R 。 
b) 寻 找 f_， 入 间 的 雅克 比 和 矩阵， 其中: 


fr 
Pe 12 
Fh Ee (12) 
Oh 
We (3) 
c) 预 测 状 态 和 矩阵 和 误差 协 方差 矩阵 
Xi = fii Xi) + XO0) (14) 
hi= FP + Qi (15) 
d) 状 态 校正 ， 计 算 卡 尔 曼 增益 矩阵 ; 


大 


HP Hi + 及 
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更 新 状态 预测 : 
X=X + KR(y — Hx i) 
估计 误差 协 方差 矩阵 : 
P=(I-KH)P, 
2.4 基于 EKF 的 速度 估计 
PMSM 的 动态 状态 方程 为 


di, 1 
TRE +PoOLi]+t 
dt LL 


Vv 
PW,Lii Po,y 1] + Re 


pi 
dt JJ2 


根据 上 述 方程 ， 可 以 写成 状态 空间 表达 的 形式 : 


B 
CA DA 


Tr di ] —R, Po,L, 0 
dt L 已 
di = w,L 一 人 -Py 4 
= - = 了 
dt L L L 
dw 多 
| | -3 -Li 3py -B 
[三 27 27 J | 
图 _ 
一 0 
L 
Vv 
二 | 0 EE | 1 
L v, 
0 0 
上 述 方程 的 离散 时 间 表 示 如 下 : 
| Po,L ] 
a oh 0 
L L 
-bo -Rh -Py 
fr Ke Ut We) = I I 7 
gq gq gq 
3P (Lb -Li 和 3Py, 7 1 B 7 
| 2J 27 " 
lr 0 
L 
Vv 
] qd 
| 二 四 
L 
0 
梯度 矩阵 为 
| Po,L 
-Er7 i 了 0 
L DL 
Pp 
FP,= fr Po,L, T 1 名 7 nr T 
6 9 L L 
-3P(L L i T 3Py, T 1 Bb T 
| 2J 2J J 


(17) 


(18) 


(19) 


(20) 


CD) 
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3 ”基于 1GSO 算法 的 EKF 参数 优化 


EKF 参数 优化 模型 

EKF 设计 的 关键 步 又 是 获得 滤波 器 参数 ， 特 别 
x(0) 、P(0) 和 协 方差 矩阵 O 、R。 算 阵 oQ 反映 了 系 
性 和 扰动 ， 和 珑 阵 尺 反映 了 由 电流 传感器 等 引入 的 疯 
本 文通 过 和 欠 代 搜索 方式 获得 这 些 系数 值 ， 从 而 获得 
计 。 改 变 协 方差 矩阵 0 和 R 将 同时 对 滤波 的 暂 态 和 
响 ; 增 大 0o 值 意味 着 


3.1 


作 期 于 
速度 估计 


是 初始 


| 


态 


统 的 不 确定 


一 


| 量 噪声 ! 


1 。 


最 优 性 能 佑 


兽 加 系统 模型 驱动 噪声 或 不 确定 性 ， 这 将 


会 加 大 状态 协 方差 项 的 值 和 滤波 器 增益 U4。 此 时 需要 提升 权重 
值 ， 使 滤波 器 的 暂 态 响应 变 得 更 快 。 类 似 的 ， 协 方差 RR 的 增加 
意味 着 有 较 强 的 衰落 噪声 , 从 而 滤波 器 的 权 值 应 该 降低 ,因此 ， 


增益 矩阵 天 的 值 将 会 降低 ， 瞬 态 性 能 变 慢 。 
对 于 初始 状态 协 方差 矩阵 P(0)， 对 角 项 表示 初 


P(0)=diagll 1 ]。 
本 文 的 主要 
试 凑 法 来 手动 选择 ， 但 这 是 一 个 非常 耗 时 的 过 程 。 


个 问题 , 本 文通 过 使 用 一 种 快速 IGSO 算法 来 调整 协 方差 矩 


始 条 件 下 的 
方差 或 均 方 误差 。 变 化 的 P(0) 产生 不 同 的 幅度 瞬 态 特性 。 但 瞬 
态 持 续 时 间 不 变 ， 且 稳 态 条 件 不 受 影 响 。 由 于 缺乏 足够 的 统计 
信息 来 评估 协 方差 矩阵 O 、R 和 P(0) 的 非 对 角 项 , 因此 假设 它 
们 是 对 角 和 矩阵 。 其 中 根据 经 验 ， 本 文 设置 x(0)=[0 0 0] ， 


标 是 选择 O 和 的 最 优 值 。 这 些 值 可 以 通过 


为 了 克服 这 


IGSO 优化 EKF 参数 的 系统 框图 如 图 4 所 示 。 
实际 状态 
Va pMSM | 
Va | 观测 输出 wn 
图 下 
| 估计 状态 ISE 
Ce ee 估计 状 
一 Q 
人 L Iso « 
J | R F = We +w,e, + We 
噪声 


图 4 ”IGSO 优化 EKF 参数 的 系统 框图 


3.2 传统 GSO 算法 


传统 GSO 算法 包含 三 个 操作 , 即 发 现 者 操作 、 搜 索 者 操作 


和 游荡 者 操作 。 在 迭代 过 程 中 ， 将 具 


最 佳 适应 度 


值 


的 成 员 选 


为 发 现 者 。 将 适应 度 值 高 于 闵 值 的 多 个 成 员 选 为 搜索 者 。 将 适 


应 度 值 低 于 阔 值 的 多 个 成 员 选 为 游荡 者 [13]。 
3.3 发 现 者 操作 


发 现 者 操作 过 程 中 ,动物 旋转 感官 受 体 从 环境 中 
在 s 维 搜索 空间 中 ， 第 z 个 搜索 回合 (迭代 ) 的 第 i 个 成 员 的 位 置 
表示 为 y; e R*， 搜索 角度 表示 为 在 = (人 ,大 ) ER ， 对 应 


和 获 信 息 


vv o 


的 搜索 方向 表示 为 F( 相 )=(f5,.…, 大 ) eR， 其 可 以 通过 极 坐 


标 变换 根据 4: 计算 得 到 ， 表 达 式 如 下 : 
A cos(%) 


fi =sin(A jn) * fi 
fi = sin(hi 1) 


( 


26) 
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彼 设 在 第 z 次 欠 代 处 的 发 现 者 位 置 为 y, ， 那 么 ， 发 现 者 将 
会 在 当前 位 置 选择 3 个 不 同 的 角度 进行 视觉 扫描 ， 即 首先 以 零 
度 扫描 ， 然 后 向 右边 扫描 ， 再 向 左边 扫描 。 设 定 视觉 的 最 大 搜 
索 角度 为 ， 视 觉 扫描 的 最 大 距离 为 4， ， 表 达 式 如 下 : 


A 7) 
i=1 


其 中 : U， 和 工 分 别 为 设计 变量 取 值 范围 的 上 界 和 下 界 。 
那么 ， 发 现 者 通过 扫描 所 发 现 的 三 个 不 同位 置 表示 如 下 : 


yzero 加 六 二 ins (4 ) 


z z 也 CO nax 
View = Vp + rida Fp (4 th (28) 


芝 党 CO nax 
Di = pt ndrax (4 1 2 


其 中 : y 表示 零度 扫描 ， ww 表示 右边 扫描 ， yw 表示 左边 

扫描 , ne R' 为 均值 为 0、 方 差 为 1 的 正 态 分 布 随 机 数 , ye Re 

为 (0,D 内 的 一 个 随机 序列 。 

然后 ， 计 算 发 现 者 搜索 到 的 3 个 新 位 置 的 适应 度 ， 并 移动 

到 具有 最 优 适 应 度 的 位 置 。 如 果 新 位 置 都 不 如 当前 位 置 ， 则 将 

其 头 转向 一 个 新 角度 ， 表 示 如 下 : 

A = A + py (29) 

其 中 : yx 表示 最 大 转向 角 。 如 果 在 a 次 迭代 结束 后 ， 发 现 者 
没有 找到 一 个 更 好 位 置 ， 则 停止 搜索 过 程 且 保持 不 动 ， 即 : 

A = G30) 


3.4 ) 搜索 者 操作 

搜索 者 操作 为 跟随 发 现 者 ,并 在 其 周围 附近 区 域 进行 搜索 。 
在 第 z 次 迭代 处 ， 第 ;个 搜索 者 根据 发 现 者 共享 的 位 置信 息 执 
行 区 域 搜索 09， 其 位 置 更 新 如 下 所 示 : 

yr = tn yi) GD 

其 中 ， eR 表示 (0,1) 区 间 内 的 随机 数 。 
3.5 ) 游荡 者 操作 

游荡 者 操作 仅 为 随机 游 走 ， 并 以 此 来 探索 新 位 置 。 如 果 将 
群 中 第 ;个 成 员 选 为 第 z 次 迭代 的 游荡 者 ， 则 它 将 生成 一 个 随 


机 角度 4 ， 其 表示 如 下 : 
二 G2) 
同样 ， 也 会 选择 一 个 随机 距离 ， 表 示 如 下 : 
d, = a nid, G33) 
然后 ， 根 据 下 式 移 向 一 个 新 位 置 : 
yr =y: +dF (4) G34) 


3.6 改进 型 GSO 算法 (IGSO) 

传统 GSO 算法 的 收敛 速度 较 慢 ， 在 对 实时 性 要 求 较 高 的 
应 用 场景 具有 局 限 性 。 为 了 使 GSO 算法 能 够 有 效应 用 到 电机 
EKF 速度 估计 系统 中 的 参数 优化 中 ， 本 文 对 其 进行 
一 种 快速 GSO 算法 ， 提 高 其 收敛 速度 。 
在 另 一 种 智能 优化 算法 , 即 粒 子 群 优化 PSO) 算 法 中 , 其 根 
据 粒子 当前 搜索 的 全 局 最 优 位 置 和 粒子 自身 的 最 优 位 置 来 更 新 


进 ， 提 出 


粒子 运动 方向 , 使 其 能 够 快速 的 向 最 终 的 全 局 最 优 位 置 靠近 5]。 
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那么 ， 可 将 这 种 思想 引入 到 GSO 中 ， 提 出 一 种 改进 型 GSO 算 
法 (IGSO), 通 过 改进 搜索 者 的 位 置 更 新 机 制 ,来 提高 收敛 速度 。 
在 第 z 次 友 代 处 ， 将 具有 最 高 适应 度 的 成 员 作为 发 现 者 ， 
将 其 位 置 标记 为 全 局 最 优 y;,，。 群 中 的 搜索 者 会 跟随 发 现 者 
在 其 领域 进行 搜索 ， 将 第 ;个 搜索 者 在 第 z 次 迭代 中 的 历史 
最 佳 位 置 标记 为 个 体 最 优 y:,,,,。 那 么 , 可 根据 全 局 最 优 和 个 体 
最 优 ， 使 搜索 者 在 跟随 发 现 者 运动 时 ， 同 时 考虑 发 现 者 位 置 和 
个 体 自身 发 现 。 这 样 就 提高 了 搜索 者 寻找 其 他 更 好 位 置 的 能 力 ， 
更 快 地 形成 发 现 者 ， 使 算法 快速 收敛 。 搜 索 者 位 置 更 新 机 制 表 
示 如 | 


1! [= 


y= thr(yopew Yt hr(yipes —y) (35) 
其 中 ，7y 为 随机 变量 ，k 和 大, 为 权重 因子 ， 这 里 都 设置 为 2。 
3.7 目标 函数 
将 EKF 模块 所 估计 的 d 轴 电流 i 、q 轴 电 流 i 和 转速 w 与 
实际 值 的 平方 误差 积分 (integral square error, ISE)， 进 行 加 权 后 


作为 目标 函数 让 ， 表 达 式 为 
F = we 十 We + We (30) 
el = [G0 六 

其 中 : ey= ,0 i) (37) 


2 
63 于 |。 加 ,1) 


正确 选择 权重 至 关 重 要 ， 否 则 可 能 会 导致 很 大 的 误差 。 本 
文通 过 进行 大 量 实验 和 结果 分 析 ， 设 置 最 佳 权 值 为 w =0.225 ， 
Ww =0.3， w=0.475。 


4 ”实验 及 分 析 


4.1 实验 设置 
利用 MATLAB/Simulink 构 建 PMSM 仿真 实验 ， 
参数 估计 算法 。 仿 真 中 的 PMSM 参数 如 表 1 所 示 。 
在 仿真 中 ， 设 置 i,i,v,v, 为 EKF 算法 的 输入 ，i,i,@, 为 
估计 状态 变量 。 为 了 模拟 实际 系统 中 的 情况 ， 在 i,,i 的 反馈 中 
加 入 了 高 斯 白 噪声 ， 强 度 为 3x10-6 ， 白 噪声 的 采样 时 间 设 置 为 
2x105 秒 。 对 于 IGSO 算法 ， 设 置 最 大 连 代 次 数 为 200， 种 群 
数量 为 20。 初始 搜索 角度 4: 取 值 为 x/4 ,最 大 搜索 角度 w， 取 
值 为 Ax/a’ ， 最 大 转向 i 取 值 为 A/2a’ ， 常 数 a 取 值 为 近似 
Vs+1 的 整数 ， 权 重 参 数 =k, =2。 


实现 EKF 


表 1 PMSM 的 仿真 参数 


参数 符号 值 

定子 电阻 ”及 0.755O 

d 轴 电 感 L 8.2mH 

q 轴 电 感 L 82mH 

磁 通 链 Wj 0.12 Wp 
转动 惯量 J 0.0025 Kg /nm 
摩擦 系数 B 0.0015 Nm/s? 

极 对 数 Db, 3 
额定 转速 @ 1000 rpm 
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4.2 参数 优化 实验 度 估计 的 精度 要 求 。 

首先 ， 进 行 验证 性 实验 ， 将 IGSO 算法 与 传统 GSO 算法 、 4.3 速度 估计 实验 
PSO 算法 和 经 典 的 遗传 算法 (GA) 进 行 比较 。 以 ISE 为 目标 函数 ， 进行 电机 启动 实验 , 设 定 启动 后 的 额定 速度 为 S00rpm。 基 


利用 各 种 优化 算法 对 EKF 的 O 和 尺 和 矩阵 进行 寻 优 ， 在 不 同 迁 ” 于 具备 最 优 O 和 RR 矩阵 的 EKF 速度 估计 器 ， 其 所 估计 的 i 和 
代 次 数 下 的 ISE 曲线 如 图 5 所 示 。 表 2 列举 了 三 种 欠 代 次 数 下 ， i 和 实际 测量 波形 如 图 6 所 示 。 由 于 状态 协 方差 矩阵 p 的 收 合 
通过 IGSO 算法 获得 的 EKF 的 go 和民 矩 阵 和 1ISE 值 。 性 问题 ， 所 估计 的 gq -g 轴 电 流 还 有 最 大 0.03 秒 的 震荡 。0.03 s 
从 图 5 可 以 看 出 ，GA 算法 的 收敛 速度 最 慢 ， 这 是 因为 其 。” 后 ， 和 矩阵 pp 收 化 此 时 i 和 i 的 轨迹 与 实际 值 接 近 。 
通过 遗传 、 交 叉 和 变异 操作 来 进行 优化 ， 过 程 较为 复杂 ， 且 容 
易 陷 入 局 部 最 优 . PSO 算法 的 收敛 速 度 比 GA 和 GSO 要 快 , 这 
是 因为 其 个 体 的 位 置 更 新 机 制 较为 优越 ， 能 够 促进 个 体 快速 寻 
找到 最 优 解 ， 这 也 是 本 文 将 该 机 制 融 入 到 GSO 中 来 构建 IJGSO 
算法 的 原因 。IGSO 的 收敛 速度 最 快 ， 大 约 在 100 次 迭代 后 就 
能 收敛 到 最 优 解 ， 而 传统 GSO 需要 约 160 次 迭代 后 才能 收敛 。 
这 是 因为 IGSO 算法 引入 了 PSO 算法 中 的 全 局 最 优 和 个 体 最 优 
融合 方法 ， 改 进 了 搜索 者 的 位 置 更 新 机 制 。 

这 就 使 IGSO 能 够 应 用 在 对 实时 性 要 求 较 高 的 优化 任务 中 。 
对 于 寻 优 所 获得 的 解 ， 当 迭代 次 数 为 100 时 ， 基 于 所 获得 的 最 
优 2O 和 丸和 矩阵 ，EKF 对 速度 估计 的 ISE 值 约 为 0.0353， 符 合 速 


[HI | 


目标 函数 值 (ISE) 


迭代 次 数 
妈 5 各 种 算法 优化 过 程 的 收敛 曲线 


表 2 IGSO 优化 的 EKF 矩阵 值 


和 迭代 次 数 对 角 怎 阵 O 对 角 拢 阵 尺 目标 函数 值 (SE) 
17.533 0 0 
1268.176 0 
50 0 6747 0 0.0451 
0 1273.562 
0 0 97.589 = 
130.637 0 0 | - ， 
75 0 6.632 0 1931:008 9 0.0375 
. 0 2334.437 
| 0 0 118.563| 
132.384 0 0 于 
100 0 7.538 0 1 0.0353 
0 2246.436 
| 0 0 126.594| 


一 一 一 ”测量 值 
一 一 估计 值 


0.1 0.2 0.3 


时 间 /s 


0.4 0.5 


图 6 启动 过 程 中 i, 和 的 估计 值 和 测量 值 


7 给 出 了 启动 过 程 中 EKF 的 速度 估计 曲线 ， 可 以 看 出 ， ”至 500 rpm, 实验 结果 如 图 8 所 示 。 可 以 看 出 , 在 额定 速度 变化 
古 计 速 度 与 实际 速度 的 匹配 时 间 只 需 0.05 s， 之 后 与 实际 速度 。 时 ，EKF 所 估计 的 速度 都 能 很 快 收敛 到 实际 值 。 
基本 一 致 。 这 说 明了 IGSO 优化 EKF 的 有 效 性 。 
为 了 进一步 验证 EKF 估计 器 对 额定 转速 变化 的 鲁 棒 性 。 设 
定 启 动 后 额定 速度 为 200 rpm， 然 后 在 0.2 s 时 从 200rpm 上 升 本 文 介绍 了 一 种 基于 EKF 的 PMSM 无 传感器 速度 控 表 


全 


5 ”结束 语 


过 
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法 , 基于 EKF 来 预测 速度 和 dq- 轴 定子 电流 。 由 于 EKF 的 性 能 
主要 取 雇 于 误差 协 方差 矩阵 O 和 R 的 值 。 为 此 ,提出 了 一 种 快 
速 的 IGSO 算法 来 优化 og 和 R ， 以 此 提高 系统 的 估计 精度 和 收 
敛 性 。 仿 真 结 果 表 明 ， 所 提出 的 方法 在 稳定 时 间 ， 噪 声 和 整体 
稳定 性 方面 表现 优异 。 

在 本 文 工 作 中 ， 电 机 负载 和 参数 都 是 一 定 的 ， 在 现实 中 ， 
环境 温度 变化 会 导致 电机 参数 变化 。 为 此 ， 在 今后 的 工作 中 ， 


将 进一步 验证 本 文 方法 对 负载 和 参数 变化 的 鲁 棒 性 。 
-1000 
0 Ol 0.2 0.3 
时 间 /s 
图 7 启动 过 程 中 转速 的 估计 值 与 实际 值 
0 0.1 0.2 0.3 0.4 0.5 
时 间 /s 

图 8 速度 变化 的 情况 下 转速 的 估计 值 与 实际 值 
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